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Abstract. We combine recent ISO observations of the vibrational ground state lines of H2 towards Photon- 
Dominated Regions (PDRs) with observations of vibrationally excited states made with ground-based telescopes 
in order to constrain the formation rate of H2 on grain surfaces under the physical conditions in the layers 
responsible for H2 emission. We use steady state PDR models in order to examine the sensitivity of different H2 
line ratios to the H2 formation rate Rf. We show that the ratio of the 0-0 S(3) to the 1-0 S(l) line increases with 
Rf but that one requires independent estimates of the radiation field incide nt upon the PDR an d the density 
in order to infer Rf from the H2 line data. We confirm the earlier result of lHabart et alJ (I2003T) that the H2 
formation rate in regions of moderate excitation such as Oph W, S140 and IC 63 is a factor of 5 times larger than 
the standard rate inferred from UV observations of diffuse clouds. On the other hand, towards regions of higher 
radiation field such as the Orion Bar and NGC 2023, we derive H2 formation rates consistent with the standard 
value. We find also a correlation between the H2 1-0 S(l) line and PAH emission suggesting that Rf scales with 
the PAH abundance. With the aim of explaining these results, we consider some empirical models of the H2 
formation process. Here we consider both formation on big (a ~0.1 fim) and small (a ~10 A) grains by either 
direct recombination from the gas phase or recombination of physisorbed H atoms with atoms in a chemisorbed 
site. We conclude that indirect chemisorption is most promising in PDRs. Moreover small grains which dominate 
the total grain surface and spend most of their time at relatively low temperatures may be the most promising 
surface for forming H2 in PDRs. 

Key words. ISM: clouds - ISM: dust, extinction - atomic processes - molecular processes - radiative transfer 

1. Introduction to this issue is to examine the H2 formation rate in differ- 

ent regions of the ISM in order to see how it depends upon 



C<3 ■ The formation of molecular hydrogen is a key process af- the local physical parameters. One can for example study 
fecting the thermal and density structure and the chem- the correlation of H 2 -related quantities (abundance, ro- 
ical evolution of the interstellar medium (see, for exam- tational excitation, vibrational excitation) with the local 
pie. ICombes fe Pineau des Foretsl l200q). Although there dust properties. ISO observations of H 2 lines and small 
is a co nsensus that H? forms on the surface of dust carbonaceous grains emission in several PDRs, and new 
grain s JGould fc Salpeterl Il963t IHollenbach fe Salpeterl H 2 UV absorption observations by FUSE of diffuse clouds, 
Il97lt lJuraL Il975l iDulev fc WilliamsL Il984|h the mecha- pen new perspectives on our understanding of the H 2 for- 
nism is not yet understood. This is partly due to our ig- ma tion process. In particular, the confrontation between 
norance concerning interstellar grain composition, form, observations and theoretical predictions p rovides stron; 
structure and physico-chemical state. It is also caused by constraints upon the H 2 formation rate faabart et al 
our lack of understanding of surface reactions in the in- |2003HGrv et all 120021: iTumlinson et~ail 12002^ . ~ 
terstellar context. j n tmg p aper) we consider PDRs and investigate the H 2 

Numerous theoretical and experimental studies have formation process using observations of H 2 emission ob- 
thus been dedicated to the study of the H 2 forma- tained by ISO and ground-based telescopes. PDRs, where 

stellar radiation plays a dominant role in determining the 
ch emical and thermal state of t he gas (for a recent review 
sec IHollenbach fc Tielensl Fl999h are privileged objects for 
the study of physical and chemical processes of the inter- 
200lt ICazaux fc Tielensl 120021 12003^1 . Another approach stellar medium. In Sect. [3 we review the data available 
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for five nearby PDRs and summarize their physical condi- 
tions. In Sect. |3 we estimate for each PDR of our sample 
the H 2 formation rate using the H 2 line intensity ratios as 
a diagnostic. Also, in order to probe the influence of PAHs 
on the H 2 formation, we compare in Sect.^the intensities 
of emission measured in H 2 fluorescent lines with those 
in the aromatic bands. Then, with the aim of explaining 
these results, we consider in Sect. [5] empirical models of 
the H 2 formation mechansim by recombination of an H 
atom in a bound site with a second H atom which is ei- 
ther in the gas phase or on a neighbouring physisorbed 
site. In Sect. \G\ we compare the observational constraints 
on the H 2 formation rate with different model predictions. 
Our conclusions are summarized in Sect.0 

2. Sample of PDRs observed by ISO 

Recently, the Sho rt Wavelength Spectrometer (SWS, 
IKessler et a,l.Lll996^l on board ISO has observed a series of 
pure H 2 rotational lines towards a variety of nearby PDRs. 
In the ISO data base, we have selected five nearby PDRs 
which sample well the range of excitation conditions cov- 
ered by the SWS observations. At the low excitation end, 
we have the Oph W and the S140 PDRs. Then, with radi- 
ation field higher by a factor of 3-10, we have the PDR IC 
63 and the reflection nebula NGC 2023 (where we focus 
on the filament at 60" south of the central star). Finally, 
at the high excitation end, we have the Orion Bar at the 
position of the peak of the fluorescent H 2 emission. These 
PDRs close to the Sun (d ~100-500 pc) are ideal targets 
to discuss the formation of H 2 in hot regions of the ISM. 

In Tabled we summarize the physical conditions pre- 
vailing in each region as determined from the literature. 
The thermal and chemical structure of the PDR depends 
on two parameters, namely, the intensity of the incident 
far-ultraviolet (FUV, 6< hv <13.6 eV) radiation fi eld and 
the ga s density. We adopt the radiation field of iDrainel 
l)l978(l and we characterize its intensity with a scaling fac- 
tor x (x— 1 corresponds to the standard FUV interstellar 
radiation field of 2.6 10~ 3 erg s _1 cm -2 ). This factor is 
determined from the expected FUV luminosity of the ex- 
citing star and assuming that the distance of the exciting 
star to the PDR is equal to the distance projected onto 
the sky. This in principle is an upper limit and hence we 
also report in Table ^ estimations of \ based on observa- 
tions of the fine structure lines of C + or 0°. In general, 
we note that x IS uncertain by a factor of about 2 to 5. 

The proton gas density uh has been derived from 
a variety of observational constraints: (1) for Oph W, 
tih ~ 10 4 cm~ 3 is estimated from the brightness profile 
of the aromatic dust emission ijrlabart et allEoO^ . while 
observations in the C° and CO lines sug gest nn ~ 6 10 4 
cm~ 3 for the inner part of the PDR ijrlabartl l200ll) : 
(2) for S140, using [C+1 158 mm and C° radio recom- 
bination line, IWvrowski et al.l l|l997b|) derive a density 
about Tin ~ 5 10 4 cm~ 3 comparable to the gas den- 
sity expected for pressure balance with the HII region 
~ 10 4 cm -3 i Timmermann et all Il99fifl : (3) measure- 



ments of CO, HCO+, HCN or CS line r atios from the 
PDR IC 63 suggest n H ~ 5 10 4 - 10 5 cm" 3 (jjansen et all 
ll994L Il995h : (4) for the bright southern emission bar 
of NGC 2023, the inferred densities are in the range 
~ 10 4 — 10 5 cm" 3 us i ng H 2 fluorescent line emission 
llBlack fc van DishoeckL Il987t IDraine fe BertoldiL Il996t 
iField et all Il998t iDraine fe Bertoldl l2000|) and ~ 10 5 
cm -3 fro m [C + 1 158 ixm and C radio re combination line 
intensity l|Wvrowski et all Il997bl l2000h : (5) finally, for 
the Orion Bar the gas density has been estimated to be 
about ~ 5 10 4 cm~ 3 fro m the observed strat ification of 
different tracers of PDRs jTielens et all Il993ft and > 10 5 
cm" 3 from obse r vation s of fine-structure line emission 
faerrman n et al lll997l), as we ll as C radio recombination 
lines llWvrowski et all Il997a| and CN, CS observations 
IjSimon et a l.. 1997). The gas density inferred from vari- 
ous atomic/molecular species shows a relatively large dis- 
persion (typically from 10 4 to 10 5 cm -3 ). This dispersion 
could result from systematic density gradients from the 
H 2 emi tting layer to the molecular cold layer (see , for ex- 
ample, IWalmslev et"afl l2000t lHabart et all l2003|) . What 
is needed for our study is the density in the H 2 emitting 
region while mainly the density tracers reflect the density 
of the cold gas in the cloud. 

We report also in Table Q the observed intensities 
for the H 2 0-0 S(3) and 1-0 S(l) lines obtained respec- 
tively with ISO-SWS and from ground based observa- 
tions. High spatial resolution observations (~1") of the 
1-0 S(l) line emission have been smoothed to the ^20" 
beam of SWS. For measurements made with beams larger 
than the SWS observations, we have scaled the flux ac- 
cording to the beam ratio. For each region, we also give 
the excitation temperature of the H 2 pure rotational lev- 
els with J < 7, the H 2 column density inferred from 
H 2 rotational lines intensity, and the 1-0 S(l)/2-l S(l) 
line ratio. In the case of NGC 2023 and the Orion Bar, 
the 1-0 S(l)/2-l S(l ) line ratio has bee n taken respec- 
tively at 18" S 11 " W lEuTton et allll99i and 50" N 30" E 
l|Walmslev et all l2000() of the SWS pointings. Finally, 
based on ISOCAM observations 1 , we give the aromatic 
dust emission smoothed to the SWS beam. In the follow- 
ing, aromatic dust particles will be hereafter referred to 
as PAHs (Polycyclic Aromatic Hydrocarbons). This is a 
generic term which encompasses large aromatic molecules 
and tiny carbonaceous dust grains containing up to a few 
1000 atoms and with radii of a few A to a few tens of A. 



3. Estimates of the H2 formation rate 

In this section, we show how the H 2 formation rate can 
be derived from the analysis of the H 2 emission line using 
PDR models. 



1 From the brightness in the LW2 filter ( 5-8.5 lira), which 
is do minated by the aromatic dust emission feoulanger et all 
ll99ST l. we can estimate the aromatic dust emission using the 
following relationship : 1 pah (2 — 15(im) ~ 2 x vl u (5 — 8.5 /im) 
based on ISOCAM-CVF spectrum (corrected from the dust 
continuum emission) taken in PDRs. 
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Table 1. Sample of PDRs observed by ISO 





Oph W 


S140 


IC 63 


NGC 2023 


Orion Bar 


x a 


250 


100 - 250 


650 


500 - 3000 


5000 - 2.4 10 4 


Ref. 


(la) 


(2a)-(2b) 


(3a) 


(4a)-(4b) 


(5a)-(5a,b) 


n H b (cm" 3 ) 


10 4 


10 4 - 5 10 4 


5 10 4 - 10 5 


10 4 - 10 5 


5 10 4 - 3 10 5 


Ref. 


(la) 


(2b)-(2a) 


(3a) 


(4c)-(4a,b,d) 


(5b,c,e)-(5e,f,g) 



A v c (mag) 10 >10 7 >10 >10 

Ref. (lb) (2c) (3a) (4a) (5 h) 

a Incident FUV radiation field expressed in units \ of the Drainc (1978) average interstellar radiation field. 
b Proton gas density nu = n H a + 2 uh 2 ■ 

c Visual extincti on within the PDR inferre d from sub-mm dus t emi ssion or CO observations . 

Refe rences: f la) IHabart et alJ EoO$); ClbllMotte et al] <1998h : (2a) IWvrowski et all (Il997rj): f2b) iTimmermann et alJ (I1996D : 
(2c) [Minchin et alJ Jl993ft: (3a) Ijansen et al] <1994 Il995l): (4 a) IWyrowski et alJ dl997bl 12000): ( 4pT TDraine fc Bertoldil J1996L 
I200Q1 : (4c) [Black fc van Dishoeckl (Il987l): (4d) iField et all Jl99 g|>: (5a) iMarconi et alJ ]l998l); (5b) iTielens fc HollenbachJ 
Jl985l): (5c) bielens et all (Il993l): (5d) iTauber et alJ Jl994h : (5e) IWvrowski et al.1 (Il997ah : (5f) iHerrmann et alJ (11997ft : (5g) 
ISimon et all (Il997l) : (5h) lHogerheiide et alJ Jl995F ! 





Oph W 


S140 


IC 63 


NGC 2023 


Orion Bar 


H 2 0-0 S(3) a 9.66 fim 
(SWS, 20") 
T^t (K) 

N C H2 (10 21 cm" 2 ) 
Ref. 


13.7[8] 

330±15 

0.6 

(1) 


16.8[30] 

500±40 

0.2 

(2) 


10 [29] 

620±45 
5 

(3a) 


16.5[13.5] 

330±15 
0.7 

(4a) 


59.7[11] 

390±20 
1 

(5a) 



H 2 1-0 S(l) a 2.12 /xm 


3.1[16] 


1.84[13] 


7[20] 


13[10] 


Beam 


1" 


74" 


1" 


1.5" 


1-0 S(l)/2-l S(l) 




2. 2 [23] 


2.8[18] 


2.3[7] 


Ref. 


(1) 


(3b) 


(4b,c) 


(5b,c) 


PAH (2-15^m)/10 3 


3.6 


0.8 


11 


62 


(ISOCAM, 6") 










Ref. 


(1) 


(3c) 


(4d) 


(5d) 



a Intensities (in 10 -5 erg s _1 cm -2 sr _1 ) with relative uncertainty in % (in between brackets). H2 line intensities have not been 
corrected for dust attenuation. The 1-0 S(l) line intensity has been smoothed to the ~20" beam of SWS or multiplied by a 
beam factor. 

b Excitation temperature of H2 pure rotational levels with J < 7. 

c H2 column density inferred from the intensity of H2 rotational lines and assuming that the population distribution of low H 2 
rotational levels is essentially in LTE. 

Refe rences: (1) Data corresponding to the peak of the H 2 emission in IHabart et all J2003I) : (2) iTimmermann et alJ dl996Tl : 
(3a) iThi et all jl999l) (note that the observ ations cannot b e fitte d by a single excitation temperature but t wo components 
at ^1 00 K and ^620 K are neede d); (3b) iLuhma n et alJ lll997t): (3c) data from Cesarsky, priv. com.; (4a) iMoutou et alJ 
lll999l): f4b)lField et all Il998ft : f4c)lBurton et all dl998f ): f4 d)lAbergel et alJ <2002h : (5a) data from Bertoldi, priv. com.; (5b) 
IvanDer Werf etHu? (119961) : (5c) IWalmslev et alJ J2000D ; f5d) ICesarskv et aITll2000D . 



For several PDRs observed by ISO, the H2 line inten- 
sities and the gas temperature as probed by the popula- 
tions of the low rotational levels of H2 w ere found to be 
higher than predic t ed by current m o dels j Bertoldi Il997 ; 
Draine fc Bertoldil Il999t IThi et all 119991 IHabart et al. . 
2003nLr^t**an ~ l2002f) . The cause of this discrepancy is 
that, in the models, the gas is not hot enough or alter- 
natively that the column density of H2 is too low in the 
zones where the gas is warm. One explanation of this dis- 
crepancy is that the H2 formation rate is larger at high gas 
temperatures, moving t he H°/H2 transition zone closer to 
the edge of the PDR. IHabart et all J2003) have shown 



that the observed H2 excitation from the moderately ex- 
cited Oph W PDR can be accounted for by increasing the 
H2 formation rate by a factor about 5 (compared to the 
standard H 2 formation rate derived in the diffuse ISM by 
|juralll975|) to ~ 2 KT 16 cm 3 s" 1 at T gas -330 K. In this 
study, we extend our study of the H2 formation in PDRs 
(described in Sect.[2J spanning a wide range of excitation 
conditions. 
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3.1. PDR model 

In order to analyse the H2 emission observations from the 
PDRs, we use an updated version of the s tationary PDR 
model described in lLe Bourlot et all l)l993|) . In this model, 
a PDR is represented by a semi-infinite plane-parallel slab 
with an incident radiation field. The input parameters are 
(i) the incident FUV field %, and (ii) the proton gas den- 
sity tih- We will consider constant density models with 
n H = 10 4 - I0 5 cm -3 and x = 10 2 - I0 5 in the range of 
values prevailing in our PDR sample (see Table QJ. With 
the inputs x an d Tin, the model solves the chemical and 
thermal balance starting from the slab edge at each Ay- 
step in the cloud. The H2 abundance results from a bal- 
ance between the formation of H2 on dust grains and the 
photodissociation of H 2 by FUV flux, which is attenuated 
by dust extinction and self-shielding in the H2 lines. At 
equilibrium, the density of atomic (n H o) and molecular 
hydrogen (nn 2 ) are given by 

R f n H n H o = R d (Q) x X e~ Td X f,(N(H 2 )) X n H2 (1) 

with Rf (cm 3 s" 1 ) the H2 formation rate, Rd(0) ~ 5 x 
10~ n s _1 the unshielded photodissociation rate per H2 
for x=l, N(H2) the colum density of H2 and e~ Td and 
f s (N(H2)) respectively the dust extinction and the H2 
self-shielding factors. In this study, we adopt Rf constant 
throughout the PDRs. The standard value inferred from 
observations of H2 UV absorption lines i n interstella r dif- 
fuse clouds is R° f ~ 3 1CT 17 cm 3 s" 1 (|.InraL Il975h . We 
assume that the energy released by the nascent molecule 
(4.5 eV) is equally distributed between the kinetic en- 
ergy of H 2 (Ek), the internal energy of H2 (Ei nt ) and 
the internal energy of the grain (E gra i n ). Moreover, we 
assume that the internal energy of the nascent H2 is dis- 
tributed in Boltzmann distribution through the energy 
levels. The energy used for the heating per H 2 formation 
is E k + E int x / with E k ~1.5 eV and E int x / -0.7 
eV where / is the fraction of the internal energy of the 
nascent H2 contributing to the heating. For the collisional 
excitation and de-excitation of H?, we a dopt the H°-H 2 
inelastic rates of lMartin fc Mandvl l)l995j) extrapolated to 
low temperatures. The heating rate due to the photoelec- 
tric ef fect on small dust grain s is derived from the formal- 
ism of lBakes fc Tielensl ( 1994). For the physical conditions 
prevailing in our PDRs, the heating is mainly due to the 
photoelectric effect. Nevertheless, for high Rf (~ 10 x R®) 
the heating rate per H2 formation becomes comparable 
(factor of ~2 lower) to the photoelectric heating rate. 

3.2. H2 line intensity ratios as diagnostic of the H2 
formation rate 

The (v, J) excited states of H2 can be populated by in- 
elastic collisions with gas phase species, UV pumping and 
by the formation process. The population distribution of 
H2 levels is a function of the gas density, the gas tempera- 
ture and the UV flux. Hence, the H2 line intensity ratios - 
which depend on the physical conditions in the photodis- 



sociation front where atomic hydrogen becomes molecular 
- should probe the H2 formation rate which controls the 
location of the H°/H2 transition zone in PDRs. In the fol- 
lowing, we study the dependence of several H2 line ratios 
as a function of Rf (predicted by model calculations). 

In Fig. ^ w e show the model predictions (face-on) for 
the intensity of the 1-0 S(l) line as well as for several 
commonly observed line intensity ratios as a function of 
the molecular hydrogen formation rate. First, we discuss 
the dependence of the 1-0 S(l) line intensity on Rf. At 
equilibrium (Eq. JU), the intensity of the 1-0 S(l) line 
Ijj is proportional to RfnjiN(H°) where N(H°) is the 
column density of atomic H atoms. For high values of 
xl n H 0.01), the H°/H2 transition is driven by the 
dust opacity and N(H°) is a constant equal to a few 
10 21 cm~ 2 . We thus have v H oc RfUu (see Fig. da)). 
Conversely, when x/ n H ^ 0.01, molecular hydrogen self 
shields sufficiently that the gas goes molecular before the 
dust gets optically thick. In this case, using the approx- 
imatio n f,(N(H 9 )) = (7V(ff2)/ An)~ 3/4 with N Q = 10 14 
cm -2 l|Draine fc Bert,oldiL ll996) for N(H 2 ) < 10 21 cm , 
we find R f n H N(H°) = 4R d (0) x x * A(ff 2 ) 1/4 A 3/4 
l|Hollenbach fc Tielensl Il999j) . At the H°/H 2 transition, 
we have N(H 2 ) = N(H)/2 with N(H) the proton column 
density from the edge of the PDR to the H°/H2 transition 
zone. Finally, oc x x A(iJ) 1 / 4 (see Fig.UJa)). 

The low-J rotational v=0-0 H2 lines have a somewhat 
different dependence upon Rf because at densities above 
10 cm -3 of interest to us, their population distribution is 
essentially in LTE and the line intensities depend mainly 
on the temperature in the photodissociation front where 
hydrogen becomes molecular. In the range of x an d n H 
studied here, the gas temperature profiles as a function of 
depth in the PDR are fairly insensitive to the fraction of 
molecular hydrogen (heating is mainly due to photoelec- 
tric emission from dust grains and cooling to fine struc- 
ture line emission of O and C + ) but higher H2 formation 
rates cause the photodissociation front to shift closer to 
the surface where the temperature is higher. We show in 
Figs. 1Kb) and[Kd) the dependence of the 0-0 S(l)/0-0 
S(3) and 0-0 S(3)/l-0 S(l) line ratios as functions of R f . 
The ratio of the two 0-0 lines is essentially a measure of 
temperature and thus decreases as Rf increases as a con- 
sequence of the temperature increase in the dissociation 
front. This effect also causes a sharper increase in 0-0 S(3) 
than in 1-0 S(l) and thus the ratio 0-0 S(3)/l-0 S(l) in- 
creases with increasing Rf. 

We also show in Fig. d c ) the dependence of the 2-1 
S(1)/T0 S(l) line predicted by the models as a function 
of Rf/RPf. This commonly observed line ratio is known 
in PDRs to vary between values of —2 typical of pure 
fluorescence to values of —3-5 at high densities (above 
10 5 cm -3 ) when collisional deexcitation of higher vibra- 
tional levels becomes c ompe t itive with radiativ e deca y 
(Bla ck fc van DishoeckL Il987t iDraine fc Bertoldil Il996j 
One sees from Fig.^c) that the model predicted line ratio 
for x < 10 3 and nn < 10 5 cm~ 3 is independent of Rf (as 
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Fig. 1. Panel (a) : Intensity of the Hi 1-0 S(l) line emission predicted by models (face-on) for two gas densities and 
different incident radiation fields - x=300 (dotted lines), x=i0 3 (solid lines) and x=l" (dashed lines) - as a fu nction 
of the value of Hi formation rate divided by corresponding to the formation rate based upon Copernicus data \,Jur<L 
1975). Panels (b), (c) and (d); several H2 line intensity ratios for the same models. 



well as of x an d nn) and equal to ^2. This is due to the 
fact that the densities and radiation fields of interest to us 
(and consequently assumed in the models) are in the range 
for which the population of excited vibrational levels of H2 
are controlled by purely radiative cascade process. 

In summary, we find that the 0-0 S(3)/l-0 S(l) line 
intensity ratio depends considerably on the H2 formation 
rate via the variation of both the gas temperature and 
the UV flux in the H°/H2 transition zone. However, the 
changes in this ratio can also be due to changes in the 
radiation field and density. Thus, by comparing this H2 
line intensity ratio predicted by the model (essentially in- 
dependent of the geometry and of the total column den- 
sity) with observations from PDRs where x an d Uh have 
been determined from other observations, we can expect 
to probe the H2 formation rate in PDRs. 

In Fig. we show the H 2 0-0 S(3)/l-0 S(l) line inten- 
sity ratio as predicted by the model for an edge-on geom- 



etry (as most of our PDRs are, see references in Tabled 
in the H2 emission region as a function of x and for differ- 
ent H2 formation rates. We compare these predicted line 
intensity ratios with observational data for our PDRs sam- 
ple described in Sect. |21 The data has been corrected for 
dust attenuation applying a correction factor r I (1 — e~ T ). 
Here, we use the extinction curve of lDrainei Jl989) and 
the visual extinction within the PDR reported in Tabled 
derived from sub-mm dust emission or CO observations. 
This assumes that our PDRs are exactly planar and edge- 
on and that densities are similar in the H2 emission layer 
and in the cold molecular layer which is questionable be- 
cause the column density inferred from the molecular hy- 
drogen rotational transitions is generally an order of mag- 
nitude lower than the value derived from sub-mm dust 
emission or CO observations. The explanation of this dif- 
ference is not clear. It could occur in an edge-on geometry 
due to beam dilution effects (i.e., H2 emission region not 
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resolved by ISO-SWS) or because of density gradients. 
For reasons of simplicity, we have adopted in Fig. [21 the 
estimate inferred from sub-mm dust emission or CO ob- 
servations. The error involved here is small since using the 
column density derived from the molecular hydrogen rota- 
tional transitions, the H2 0-0 S(3)/l-0 S(l) line intensity 
ratio would diminish by ^30% which is comparable to the 
error bars. 

In the following, we discuss for each PDR of our 
sample the value of Rf required to account for the ob- 
served H 2 excitation. For Oph W (x ~250, n H ~10 4 
cm" 3 ) and IC 63 (x -650, n H ~ 5 10 4 - 10 5 cm" 3 ), 
we find that models with a high H2 formation rate (i.e., 
R f > 5 x R° f ) roughly reproduce the 0-0 S(3)/l-0 S(l) 
line ratio observed. For the standard H 2 formation rate, 
the 0-0 S(3)/l-0 S(l) line ratio is underestimated by a fac- 
tor about 4. For the more highly excited PDRs, i.e., NGC 
2023 (x ~ 500 - 3000, n H ~ 10 4 - 10 5 cm" 3 ) and the 
Orion Bar (x ~ 0.5 - 2.5 10 4 , n H - 5 10 4 - 3 10 5 cm" 3 ), 
we find on the contrary that models with the standard 
H2 formation rate roughly match the data. In the case 
of S140, where we have no measurement of the intensity 
for the 1-0 S( l) line, we use ISO obser vations of the 1-0 
Q(3) line from lDraine fc Bertoldil {l999) to determine the 
H2 rotational to rovibrational line intensity ratio. For the 
physical conditions of interest to us, the models predict 
that the 1-0 S(l)/l-0 Q(3) line intensity ratio is about a 
factor of 1.4. Applying this factor, we find that the 0-0 
S(3)/T-0 S(l) line intensity ratio (corrected for dust at- 
tenuation) is about 3. Then, one sees from Fig.[5]that for 
S140 (x ~100-250, n H ~ 10 4 -5 10 4 cm" 3 ) models with a 
high H2 formation rate are - as for the moderately excited 
PDR Oph W and IC 63 - required to explain the observed 
H2 line intensity ratio. In Table |3J we give for each region 
the values of the H2 formation rate derived using Fig. 
and assuming that the gas density in the H2 emission re- 
gion is about 10 4 cm~ 3 for all sources. Considering that 
the gas density is about 10 5 cm -3 , Rf required to explain 
the data would be (by a factor of ~2) higher and lower for 
respectively the moderately and the highly excited PDRs. 
We emphasize that our determination of Rf from the H 2 
line intensity ratios depends on the treatment in the model 
of the thermal balances and the chemistry. In particular, 
the calculation of the grain photoelectric heating which 
mainly determines T gas in the PDR modelling (see Sect. 
13. 1J) could affect the excitation of the lowest rotation al lev- 
els (J <5) of H 2 . For example. lHabart et alJ l)2003h have 
shown that, throughout the Oph W PDR, increasing the 
photoelectric heating rate by ^50%, which leads to higher 
T gas in the H°/H2 transition, the 0-0 S(3) line intensity 
is enhanced (by a factor ~2) and Rf required to explain 
the data would be reduced by a factor ~ 1.5 (see Fig- - 
However, even if we cannot precisely determine the uncer- 
tainties on our estimates of Rf, the values of Rf derived 
here should be at most uncertain by a factor of about 2. 

One check of our estimates of Rf is to see if the pre- 
dicted T gas values in the H2 emitting region are consis- 
tent with the observed rotational temperatures reported 



in Tabic [J In Fig. EJ we show the 0-0 S(3)/l-0 S(l) line 
intensity ratio predicted by the same models shown in Fig. 
El as a function of the gas temperature in the H2 emitting 
region and compared with observational data. For obser- 
vations, we use the excitation temperature of H2 pure ro- 
tational levels (with J < 7) which in principle is an upper 
limit as UV pumping could contribute to the excitation of 
H2 even for low energy levels. In the case of IC 63 where 
the measured pure H2 rotational lines sh ow two exc i tation 
temperatures (at ~100 K and ~600 K, iThi et all Il999j) 
there is direct evidence for such contamination. Using Fig. 
13 we favour (as previously) for the moderately excited 
PDRs models with high H2 formation rate which can in 
fact explain the observed rotation temperatures. Models 
with a standard H2 formation rate predict in fact a gas 
temperature (<300 K) lower than observed (^300-600 K). 
For the Orion Bar and NGC 2023, we find on the contrary 
that models with standard H2 formation rate predict tem- 
peratures consistent with the data (~300-400 K). Finally, 
we have checked that these models reproduce the observed 
absolute intensities of the 0-0 S(3) and 1-0 S(l) H2 lines. 
Taking into account the inclinations of the PDRs of our 
sample, which are seen edge-on except for IC 63, we find 
an agreement within <50%. 

It must be emphasized that these results are based on 
two fundamental assumptions. Firstly, we assume a uni- 
form homogeneous gas density while PDRs may have a 
density gradient and could be clumpy. Nevertheless, be- 
cause the modelling of the H2 excitation mainly depends 
upon the average gas density in the H2 emission zone, 
density structur e effects are probabl y minor. In fact, the 
detailed study of lHabart et all l|2003j) of the H 2 excitation 
toward Oph W taking into account the gas density profile 
reaches the same estimation of Rf deduced here. 

Secondly, we assumed a static, equilibrium PDR. In 
reality, the propagation of the ionization and photodis- 
sociation fronts will bring fresh H2 into the zone emit- 
ting line radiation. Note that a non-equilibrium ortho- 
to-para H2 ratio has been observed in I SO-SW S ob- 
servat ions towards Oph W l|Habart et all 12003ft. NGC 
2023 l)Moutou et allll999|) and NGC 7023 l)Fuente et all 
199 91) . Nevertheless, the mode l of non-equilibrium PDRs 
o flStorzer Hollenha,chl |l998) predicts that, for the phys- 
ical conditions prevailing in the Orion Bar and for an 
advection velocity of the order of 1 km s" 1 , the H2 0-0 
S(3)/T0 S(l) line intensity ratio varies by a factor less 
than ~2 relative to the steady state value. 

In the next section, we review other observational con- 
straints on the H2 formation rate based on H2 UV absorp- 
tion measurements. 



3.3. Other observational constraints on the H2 
formation rate 

UV absorption measurements of H2 have been often used 
to study H2 formation, destruction, and excitation in 
the diffuse ISM. With the numerous new H 2 UV ab- 
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Fig. 2. H 2 0-0 S(3)/l-0 S(l) line intensity ratio, as predicted by models (edge-on) and as observed for our PDR sample 
(crosses with error bars), as a function of x the FUV incident radiation field. Predicted line intensity ratios (taken 
in the H2 emitting zone) are presented for two gas densities and shown for three different H2 formation rates (R°f, 
5 x BP., 10 x Rfj). Observed line intensity ratios have been corrected for dust attenuation (see text). For each PDR of 
our sample, x has been derived from the expected FUV luminosity of the exciting star and assuming that the distance 
of the exciting star to the PDR is equal to the distance projected onto the sky (see Sect. Wjl- For S140, where we have 
no measurement of the intensity for the 1-0 S(l) line, we use observations of the 1-0 Q(3) line to determine the H% 
line intensity ratio (see text). 



sorption line obser v ations obta i ned re c ently by F USE 
1 1 1 1 1 1 1 1 1 1 1 1 1 



lis now et all 


2000; 


Shnll et all l2000t iRachford et al. 


2001, 
20021) 


2002; 


Tumlinson et all 2002: Sonnentrucker et al. 


the H 2 


formation in the diffuse ISM has been re- 



considered. In particular, iGrv et a <|2002l) have deter- 
mined the H2 formation rates over three lines of sight 
in the Chamaeleon clouds. They find a rate roughly con- 
stant and equal to about 4 10 -17 cm 3 s _1 (with an un- 
certainty of about a factor of 2) for riu ^ 30 — 50 cm~ 3 



an d T nn il ~ 60 K, in agreement with the rate inferred 
bv l.lural <)l975h . Moreover, due to the high FUSE sensi- 
tivity, fainter stars with higher extinctions could b e ob- 
served with far-UV instruments l|Moos et al||2000 |). thus 
allowi ng the study of translucent clouds. IRachford et all 
ll2002l) . studying correlations of the HVrelated quantities 
with the column densities of other molecules and dust 
extinction properties for lines of sight with A v > 1, in- 
vestigate the dependence of the H2 formation rate with 
the composition and physical state of the gas and grains. 
Furthermore, the FUSE observations in the lower metal- 
licity environments of the Small and Large Magellanic 
Clouds allow us to probe H2 formation and destruction 
in phys ical and chemical enviro nments different from the 
Galaxy. iTumlinson et alJ (J2002) find that to reproduce the 
reduced molecular fraction and enhanced rotational exci- 



tation in the SMC and LMC, a low H2 formation rate 
(Rf ~ 3 10~ 18 cm 3 s^ 1 ) and a high UV field relative to 
diffuse Galactic medium are required. 

Combining the FUSE results in the diffuse 
Chamaeleon clouds and the ISO observations dis- 
cussed above, we find that for a wide range of physical 
conditions - l<x< 10 4 , 100 cm" 3 < n H < 10 5 cm" 3 , 
50 K< T gas <600 K, 10 K< T dust <100 K (see Tables Q] 
and [3J| - H 2 forms efficiently {Rf -4 10" 17 - 1.5 10" 16 cm 3 
s _1 ). This result raises questions about our understanding 
of the H2 formation process. In Sect. [SJ we re-examine 
models of H2 formation with this in mind and compare in 
Sect. El the observational constraints on the H2 formation 
rate with model predictions. 

4. Influence of aromatic dust on the H2 formation 
process in PDRs 

In PDRs, small dust grains are intimately coupled to 
the evolution of the gas. In fact, recent theore tical 
llBakes fc Tielenslll994llWeingartner fc Drainell200l|) and 
observational l|Habart et all l200l|) work has shown that 
small grains (radius < 100 A) dominate the photoelectric 
heating. Furthermore, given that small grains make a dom- 
inant contribution to the total grain surface (see Sect.EJ), 
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Fig. 3. H 2 0-0 S(3)/l-0 S(l) line intensity ratio as predicted by the same models shown in Fig. and as observed for 
our PDR sample (crosses with error bars) as a function ofT gas the gas temperature in the H2 emitting region. For the 
observations, we use the excitation temperature of the H2 pure rotational levels (see text). For IC 63, the observed H2 
pure rotational lines intensity cannot be fitted by a single excitation temperature but two components at ^100 K and 
^600 K are needed. 



it is plausible that they play a dominant role in H2 forma- 
tion. Therefore, if H2 forms on small grains, we can expect 
that both the grain photoelectric heating rate and the H2 
formation rate will scale with their abundance. In partic- 
ular, increasing their abundance would lead to higher gas 
temperature and would bring the H°/H2 transition closer 
to the edge: this could significa ntly enhance the rota tional 
and fluorescent emission of H2 . lHabart et alJ l)2003[) have 
shown that, assuming that Rf scales with the PAH abun- 
dance, the enhancement in Rf required to account for the 
observed H2 emission from Oph W may result from an 
increased PAH abundance. Based on this, we investigate 
what we can predict concerning the influence of PAHs on 
the H 2 formation process. 

To free oneself of the influence of the PAHs on the 
thermal balance, we study here the ratio between the H2 
1-0 S(l) line (whose intensity does not depend on T gas , see 
Sect. I3.2)l and the PAH emission. We first express analyt- 
ically the H2 fluorescent to PAH emission ratio I H / 1 pah 
as a function of the H2 formation rate and the PAH abun- 
dance, as well as the physical conditions x), in order 
to examine qualitatively the dependence of this ratio with 
these parameters. In Sect. 13.21 we have seen that the H 2 
1-0 S(l) line emission Ijj goes roughly as (i) RfUn for 
high x/ n H and as (ii) \ f° r l° w x/ n H- The emission of 
aromatic dust scales with the intensity of the FU V radia- 
tion field l|Puget et allll985HSelkren et allll985|) and we 
have Ipah oc x e" Td x N H [C/H\ PAH with [C/H] PAH the 



abundance of carbon locked up in PAHs. As the column 
density over which PAHs emits is a few 10 21 cm" 2 (where 
dust is optically thin), Ipah goes as % x [C/H]pah- From 
these considerations, we find 



L H 2 



LPAH 



[C/H]_ 



II H 

X 



[C/H]PAH 



X /n H > 0.01 
X /n H < 0.01 



(2) 



Thus, by studying the observed H2 fluorescent to PAH 
emission ratio from PDRs covering a wide range of nn/x 
ratio, we should be able to probe the H2 formation on 
PAH: if H 2 forms on PAHs, we expect that the I H2 /Ipah 
ratio scales with n^/x- In other words, the Rf /[C / H]pah 
ratio is constant and does not vary from one region to 
another. 

We now compare the ratio between the H2 1-0 S(l) line 
and PAH emission observed from the PDRs of our sample 
to model results. For the models, we consider x = 10 2 — 10 4 
and nn = 10 4 — 10 5 cm~ 3 which corresponds to the range 
of values prevailing in our PDR sample (see Table We 
adopt a constant Rf /[C/ H]pah ratio corresponding to 
the values derived in the Oph W PDR, i.e. , Rf ~ 5 X 
R° f and [C/H} PA h -0.5 10" 4 llHabart et all 1200.1) . The 



power emitted by PAHs h as been derived from the PAH 
absorption cross-section of lVerstraete fc Leeerl l|l992i) and 
the size distribution described in Sect. In Fig. we 
show the H2 1-0 S(1)/PAH emission ratio predicted by the 
model in the H2 emission zone as a function of the nn/x 
ratio. As expected from the Eq. @, for x/ n H ^ 0.01 the 
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0.0100 



I 

< 



0.0010 




0.0001 



Fig. 4. Ratio between the H2 1-0 S(l) line and PAH emis- 
sion as predicted by edge-on model (solid line) and as ob- 
served (crosses with error bars) as a function of nn/x- 
The dotted straight lines shows the linear dependence of 
the emission ratio with njj / \ expected from Eq. 0) when 
x/ n H ^ 0.01. For the PDRs of our sample, nn/x has 
been taken equal to the ratio between the average value of 
the uh determinations given in Tableland \ determined 
from the expected FUV luminosity of the exciting star and 
assuming that the distance of the exciting star to the PDR 
is equal to the distance projected onto the sky (see Sect. 



formation by physisorbed H atoms (Langmuir- 
Hiiishclwood) whereby two adsorbed and mobile ph- 
ysisorbed H atoms (bound simply via van der Waals 
interaction w ith binding energi es of the order of 0.05 
cV or 500 K, iKatz et all Il999h interact to recombine 
and desorb as H2 and 

formation by the interaction of an H atom from the gas 
phase with a chemisorbed H atom (wit h binding energy 

10,000 K. lFromherz et all 



roughly 1 eV equivalent to ~ 
1993) forming desorbed H2 
nism) . 



(the Eley-Rideal mecha- 



In the following discussion, we describe these pro- 
cesses mainly focusing on the version of the Langmuir- 
Hinshclwood mechanism whereby one of the ad- 
sorbed H-atoms is or i ginally in a chemisorbed site 
llHollenbach & SalneterL Il97ll ICazaux fc Tielensi 120021 
2003). This approach will be called in the following the 
indirect chemisorption approach as opposed to the direct 
chemisorption or Eley-Rideal mechanism. We make sim- 
ple empirical estimates of rates for these processes as a 
function of the physical conditions and grain characteris- 
tics. We consider small grains (SGs), dust particles with 
radii of a few A to 100 A (i.e., PAHs and Very Small 
Grains, VSGs) which fluctuate in temperature in the ra- 
diation field, and big grains (BGs), larger grains of radii 
>0.01 urn in t herm al equilibrium with the radiation field 
ilDeserl el all fl<)<)il,. 



H2 1-0 S(1)/PAH emission ratio increases proportionally 
with nn/x- 

The observational ratios for our PDR sample are com- 
pared to these predictions. The observed points roughly 
fall on the model curve. From one PDR to another, the 
change of the emission ratio seems to result mainly from 
the physical conditions (i.e., nn/x) variations. We deduce 
from this that the Rf /[C/ H]pah ratio is roughly con- 
stant which suggests that formation of H2 on PAHs should 
be important. However, the correlation between Rf and 
[C / H]pah found in our PDR sample is observed not to 
apply to the main PDR in the 30 Dor star forming re- 
gion in the Large Magellanic Cloud where for nn/x ~ 1 
the H2 1-0 S(l) line to PAH e missio n ratio is measured to 
be ~0.001 (|Boulanger et all l2003t) . In future work, this 
comparison should be extended so as to include all small 
grains, i.e., not only the band carriers but also the very 
small grains emitting at longer wavelengths. 

5. Ho formation mechanisms on grain surfaces 

In this section, with the aim of explaining our observa- 
tional constraints on the H2 formation, we examine mod- 
els of H2 formation mechanism. Two general mechanisms 
for forming H2 on grain surfaces have generally been pro- 
posed: 



5.1. The Langmuir-Hinshelwood mechanism 

In the formation from physisorbed atoms (pure Langmuir- 
Hinshelwood), an H atom is already physisorbed on the 
grain surface. A second H atom from the gas sticks to 
the grain and diffuses to finally recombine with the first 
H atom to form H2. The rate depends upon the compe- 
tition between the mobility of the H atom on the grain 
surface and its thermal evaporation rate. In the labora- 
tory, this process is observed to be efficie nt on g rains with 
7 < T dust < 20 K JPirronello et a ll Il997l ll 99Ql) . 

Amongst the dust grain populations described above, 
only BGs could carry a large fraction of physisorbed H 
atoms. Indeed, after a UV photon is absorbed by SGs, 
the physisorbed H atoms should evaporate. Further, for 
the physical conditions typical of the PDRs, the rate of 
thermal fluctuations t~^ s (~ 1.6 10~ 9 x N(C) x x s -1 
with N(C) the numb er of carbon atom in the grain, 
IVerstraete et al 1 12001 is comparable to or larger than 
the accretion rate of H atoms t^ c : for a grain with a 
radius of 1.5 nm (N(C) ~ 1600), we find t'^Jt^. ~ 
1000 x (100 K/T gas ) - 5 x x/n H - Thus for X /n H larger 
than 0.001 (true in general for the objects in Table rj), 
UV radiation in PDRs will keep SG surfaces clean. 

Since this mechanism will only be efficient on BGs and 
in interstellar regions of low excitation (in PDRs we find 
typically Tbgs > 30K, see Table EJ, we neglect it in what 
follows. 
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5.2. Formation involving a chemisorbed H atom 

We now consider the case where the first H atom 
is bound to the surface in a chemisorbed site. 
The second H atom from the gas phase reaches 
this site dire ctly (direct chemisorp t ion o r Eley-Rideal, 



Table 2. Parameters for carbonaceous dust surface. 



iDulevt Il99& IParneix fe Brechignad Il998|) or after dif- 
fusion on the grain surfac e (indirect chemiso r ption, 
Hollenha, ch fc Salneieii Il97lt ICazanx fc Tielensi 120021 
2003). The incident H atom crosses the activation bar- 
rier E„ (typically about ~50 -1500 K. iFromherz et al 



19931 IParneix fc Brechignad ll99St ISidis et all I200C 
Cazaux fc Tielensi 120031) by either thermal hopping or 
tunneling and recombines to finally form H2. This pro- 
cess involving strongly bound H atoms will be efficient at 
higher dust surface temperatures than the pure Langmuir- 
Hinshelwood mechanism. 

We first consider the direct Eley-Rideal mechanism. 
This involves reaction of the incoming H atom with an 
chemisorbed H on the grain surface without any require- 
ment for adsorption and thermal accommodation. The for- 
mation probability Pjj 2 can be written as: 

Ph 2 = f x ry (3) 

with / the probability that the incident H atom reaches 
the chemisorbed H and recombines to form H2 and r\ the 
probability to eject the molecular hydrogen formed. The 
incident H atom from the gas phase must hit the grain 
with enough thermal energy to cross the activation barrier 
(E a ) by thermal hopping. / is given by : 

E a 



t = — x exp 
J N y 



kTgas 



(4) 



with N c the total number of chemisorbed H atoms on the 
grain surface and N the total number of physisorbed and 
chemisorbed sites. The fraction of chemisorbed H atoms 
(N c /N) is critical for this process. This fraction will be 
significant when T gas is sufficiently high to carry an ap- 
preciable fraction of chemisorbed sites occupied. In PDRs 
(where T gas >300 K) we expect N c /N to be relatively 
high while in cold interstellar clouds N c /N should be low. 

An alternative to the above is the case where the inci- 
dent H atom initially sticks to a physisorbed site and then 
scans the surface until it "finds" the chemisorbed H atom. 
Here, the formation probability Ph 2 can be written as: 

Ph 2 = S x / x r, (5) 

with S the sticking probability. 

The probability / that the physisorbed H atom diffuses 
over the grain surface (until it reaches a neighbouring site 
to a chemisorbed H atom) and recombines before it is 



evaporated, can be given by / = 



where r_ is the 



inverse time scale for diffusing and recombining and 
the evaporation rate of physisorbed H. In the appendix 
we estimate the time scales t v and r ev . This leads to the 
expression: 

1 



/ = 



1 + (^)* x exp + ( fa) x k r 



(6) 



Parameter 


Value 


Ref. 


E d (K) 


600 


(1) 


E P (K) 


500 


(1) 


E a (K) 


600-1000 


(2)-(3,4) 


Ax (A) 


~2 


(2) 



Ed is the desorption energy of a physisorbed H atom; E p is the 
activation barrier energy for physisorbed H diffusion; E a and 
Aa; are respectively the energy and the width of the activation 
barrier for recombination with a neighbouring chemisorbed 
H. Here, we assume that the barrier E a for recombination to 
form H2 is the same between a physisorbed and a chemisorbed 
site as for direct recombinati o n from the ga s phase- 
Refere nces (1) iKatzetalJ (Il999lh (2) ICazaux fc Tielens 
1200.^: (■{) IFromherz et alJ ijlflfllb f4l IParneix fc Brechignac 
(1998) 



with Ed the desorption energy of a physisorbed H atom 
and E p the activation barrier energy for the diffusion of 
a physisorbed H (see Tabled- Considering that the ph- 
ysisorbed H atom crosses the activation barrier to recom- 
bine with a neighbouring chemisorbed atom by either (i) 
thermal diffusion or by (ii) tunnelling, we find (see ap- 
pendix) that k rec — exp \ a t high dust temper- 
atures (above ~40 K) when thermal diffusion dominates 
whereas at low temperatures (below ~ 30 K) when tun- 
neling is more important k rec — exp PE^(jgT- 1 



kT du 



with kT cr = 80 E d3 /(E°- 3 5 x Ax(A)) where E a3 = 
£ a /(1000 K), E d3 = £ d /(1000 K) and Aa; is the width of 
the barrier typically of the order of ~l- 3 A llBuchl Il989t 
IFromherz etail Il993t ISidis et alJ.l2000j) . 

The activation barrier energies for physisorbed H atom 
diffusion (E p ) and for recombination with a neighbouring 
chemisorbed atom (E a ) are critical for this process. For 
E a high (>1000 K), the recombination term will dominate 
and / will be high only where tunneling dominates (see 
Fig. 0. For E a low (^600 K), the diffusion term will be 
critical and / will depend much on E p . In this case, / will 
be high for E p < Ed (if E p is too large the atoms evaporate 
before they find a neighbouring chemisorbed site) and for 
E p ~ 500 K we find that the fraction of chemisorbed H 
(N c /N) needs to be at least about 0.1 (see Sect.EJ). 

We note that formation by chemisorption can hap- 
pen on the surface of big or small grains. In fact, the 
thermal fluctuations undergone by SGs after absorbing 
a UV photon are not sufficient to evaporate chemisorbed 
H atoms (binding energy around 10,000 K). Moreover, we 
estimate that the timescale for a physisorbed H to find a 
chemisorbed H is much less than the timescale of thermal 
fluctuations (see Sect. 15. l|) . 

Since SGs make a dominant contribution to the to- 
tal grain surface (see Sect. EJ) and also because they may 
have numerous chemically bonded hydrogen atoms (be- 
cause of their small size they are more disordered), it is 
plausible that they play a dominant role in H2 formation 
by chemisorption. However, if small grains contribute effi- 
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ciently to H2 formation, they must be continuously rehy- 
drogenated : this could be the case for larger PAHs and 
VSGs for which the adsorbed H atoms should evaporate 
less after UV photon absorptio n because of the lower te m- 
perature fluctuations (see, e.g.. lVerstraete et all l200lj) . 

In the next section, we compare the H2 formation rates 
as predicted by the different mechanisms discussed here 
and as derived from observations. 

6. Comparison between the H2 formation rate as 
predicted and as derived from observations 

In this section, we see how the observational values of H2 
formation rate in PDRs fit and constrain H2 formation 
mechanisms on grain surfaces. 

First, based on the empirical formation probabilities 
given above, we estimate the rate of formation through 
chemisorption on small and big carbonaceous grains. The 
formation rate Rf (in cm 3 s _1 ) can be written as: 

R f = ^av H x P„ 2 (7) 

where avjj is the collision frequency between H atoms 
(with mean velocity Vh) and grains and Ph 2 the proba- 
bility of formation (see Eqs. and JSJ). The mean cross 
section for collisions between grains and H atoms, a (cm 2 ), 
is given by < " a ™'^ ^ ^a 2 > with n gra i n (a) the grain den- 
sity with radius between a and a+da. We assume a MRN 
size distribution with a=3.5 while the lower and upper 
limits of the grain radius are respectively a m i n =A and 
a m ax=W0 A for SGs and a m in=0-0l and a max ^0.1 fim 
for BGs. Adopting uniform mass density p=2.25 g cm~ 3 
(typical value for graphite grains) and an grain/gas mass 
ratio G about 0.001 2 for SGs and 0.01 for BGs, we esti- 
mate 3 : (i) for SGs, av H ~ 6 10" 16 (^f^) a5 cm 3 s" 1 ; (ii) 
for BGs, av H ~ 3 10~ 16 ( ^fn< ) a5 cm 3 s" 1 . Under these 
assumptions, we find that the collisions cross section is 
—2 higher for SGs than for BGs and that there is enough 
grain area with SGs to make Rf larger than the standard 
H2 formation rate. We assume the small dust abundance 
to be constant from one region to another although vari- 
ations are expected. However, this assumption should not 
be critical for our approach which attempts with simple 
empirical models to see what type of H 2 formation pro- 
cesses are to a first order relevant in PDRs. 

To calculate Ph 2 j we adopt the following assump- 
tions. Firstly, for the indirect chemisorption, the for- 
mation efficiency will depend on the sticking coeffi- 
cient S which remains uncertain. Studies of the stick- 
ing of H on grain surfaces generally predicts that S 

2 The carbon locked up in SGs has an abundance of 
[C/H] ~ 10 -4 inferred from comparison between observations 
of du st galactic emission and extinction wit h model calcula- 
tions JPesert et all Il990l: lLi~&: Drainell200lfl . 



is about 1 at low temperatures and decreases with in- 



We calculate ovh with vh 



,SkT, 

\ TT 



^ a min 



"da x 



with 



vl/2 



and 



faollenbach & McKeel lie 


7& 


p„ + "-^ ^ j 

Burke & Hollcnbach, 


llQSsl 


iLeitch-Devlin & Williams 


1198,4 iBuch & Zhand. 


1991D 



V, 



grain / gas 



fa max ~5' ITa ^ a a da 

1.4 tuh G/p the volume of grain per hydrogen nuclei 



Here, we approxi mate S by 1/(1+ T nns , / 400 K + (T nnst 
I 400 K) 2 ) from iBurke fc Hollenbacbl Jl983lh iBertoldil 
Jl997h . Secondly, for the parameters of the physisorbed 
and chemisorbed sites, we assume typical values for car- 
bonaceous grains which are summarized in Table [3 We 
treat the surface coverage of chemisorbed H atoms (N c /N) 
as a free parameter. Further, we take N c /N for SGs and 
BGs to be equal and constant although this fraction is 
likely to depend upon grain characteristics and physical 
conditions. This is not completely satisfactory but it gives 
some insight into this parameter which more detailed the- 
ories of these processes should attempt to fit. We find that 
N c /N ~0.1 (see below) gives a reasonable fit to our PDR 
data. Finally, we take the probability to eject H2, 77, to 
be unity. This is based simply on the fact that ejecting 
the newly formed H2 molecule (the desorption energ y is 
-0.05 eV, iKatz et all ll QQflt iGazaux fc Tielenst 1200,^ re- 
quires a small fraction of the available formation energy 
deposited in the grain (1.5 eV). Note that if we do not as- 
sume spontaneous desorption, the H2 formed would stay 
on the surface for really low temperatures (below 10 K) 
since the t emperature is not high enough to allow evapo- 
ration (see ICazaux fc Tielensl [20021 120031) . 

For the purpose of comparison with the formation rates 
determined using PDR models (Sect.|3J), we need to esti- 
mate dust temperatures in the PDRs of interest to us. 
In Table |3 we give the temperatures of SGs and BGs 
expected in the H°/H2 transition zone, where the radia- 
tion field is given by \ ex P ( — T ) with r the FUV dust 
opacity in the transition zone equal to about unity. For 
the big grains at thermal equilibrium, we h a ve ad opted 
the analytic expression of iHollenbach et al ] l)l99llh For 
the small grains subject to thermal fluctuations, we list 
in Table the temperatures for a graphite grain of ra- 
dius 1.5 nm (N(C) ~ 1600). We have computed for each 
value of the radiation field the full temperature histogram 
and the values in Table correspond to the median value 
of the time dependent temperature (i.e., the grain spends 
half of the time at temperature higher and lower than 
this value). The median temperature of small grains is 
unlike that of big grains not constrained by their emis- 
sion spectrum. In these calculations, we used the absorp- 
tion and emission cross section of bulk graphite from 
Draine fc Led (119841) : the heat cap a city is also of graphite 
given by iGuhathakurta fc Draind <|1989j) . We have thus 
ignored quantum effects which might significantly change 
properties of these particles. This estimate of the small 
grain temperature is only indicative. Ideally, one should 
take into account the time dependence of the tempera- 
ture as a function of the grain size. Further, their emis- 
sion properties might not allow then to cool down to the 
cosmic 2.7 K background temperature. Moreover, for very 
low temperatures (< 10K), the heating due to the colli- 
sions between the grains and warm gas {T gas >100 K) or 
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Table 3. Temperatures of dust and gas and H2 formation 
rates 



Region 


rpa 

J SGs 


1 BGs 


rpb 

1 gas 


Rf c 




(K) 


(K) 


(K) 


(cm 3 s" 1 ) 


Chamaeleon 


>2.7 


15 


60 


4 10" 1 ' 


Oph W 


10 


36 


330 


1.5 nr 16 


S140 


10 


36 


500 


1.5 10~ 16 


IC 63 


12 


44 


620 


1.5 10~ 16 


NGC 2023 


25 


60 


330 


3 10 -17 


Orion Bar 


62 


90 


390 


3 10 -17 



a Temperature of the small grains (SGs) and big grains (BGs) 
expected in the H°/H2 transition zone (see text in Sect. [(J. 

Gas temperature inferred from the distribution of low H2 
pure rotat i onal levels (see Table 0. For the Chamaeleon see 
iGrv et al.l d2002l) . 

c H2 formation rates derived from observations (see Sect.^Jl. 



due to the formation of H2 (—1/3 of the formation energy 
is deposited in the grain) would increase the grain tem- 
perature. However, in the following we will use the SG 
temperatures reported in Table El in order to see if small 
grains are relevant for forming H2. We find that the me- 
dian temperatures of small grains are significantly lower 

than those of big grains. Indeed, SGs spend most of their 

1 1 1 1 

time at low temperature ( Guhatha kurta & DraineL ll989'). 

With the above formalism, we compute the formation 
rate Rf for the various processes outlined above on both 
SGs and BGs. The results are presented in Fig. El where we 
compare the H2 formation rate predicted for indirect and 
direct chemisorption with the rates derived from observa- 
tions. For completeness, we also show the results for the 
case of the pure Langmu ir-Hinshelwood me chanism upon 
BGs using the model of iKatz et al.1 1 1999h based on ex- 
perim ental results for graphite surfaces ijPirronello et all 
119991) . 

First, we see from Fig. E{a) that indirect chemisorp- 
tion may explain the data for all sources if the activation 
barrier energy for recombination is E a < 1000 K and if 
the fraction of chemisorbed H atoms is as large we have 
assumed, N c /N ~0.1. This result is broadly consistent 
with t he re cent theoretical model of ICazaux k. Tielensl 
(2002, 2003), which accounts properly for the population 
of chemisorbed and physisorbed sites on the surface. Thus, 
our estimate of the fraction of chemisorbed H atoms may 
be reasonable. However, the value of the activation barrier 
energy is uncertain and critically influences our estimate 
for Rf and the temperature range for which formation 
via indirect chemisorption is important. We find in fact 
that for E a — 1000 K this process will be efficient until 
Tdust ~30 K (i.e., important in PDRs of relatively moder- 
ate excitation, x <5000); for Tdust >30 K the physisorbed 
H atoms evaporate before recombining. Conversely, for 
E a —600 K the formation via indirect chemisorption will 
be still efficient at Tdust —100 K (i.e., important even in 
highly excited PDRs). In this last case, where the diffusion 
term dominates (see Sect. I5.2jl . the decreasing of Rf from 



moderate to highly excited PDRs is due to the competi- 
tion between evaporation and finding the chemisorbed H 
atoms over the surface. 

Secondly, we find that small grains which dominate 
the total grain area and spend most of their time at low 
temperature (< 30 K for x <3000, see Table EJ) may be 
the most promising surface for forming H2 via indirect 
chemisorption. The H2 formation on SGs can be effective 
even in highly excited PDRs where the equilibrium tem- 
perature of BGs can be too high to form H2. In the case 
of SGs, we note that in most objects, the observed H2 
formation rates can be explained with a high value of E a 
(-1000 K). For BGs, conversely a low value of E a (<600 
K) is required. However, our description of the temper- 
ature of SGs is crude (median value of the temperature 
distribution for a single size of grain). As a consequence, 
we cannot quantitatively constrain the parameters (E a , 
E p , ..) of the H2 formation process. 

The direct or Eley-Rideal mechanism fails to explain 
the observations by a factor of a few (see Fig. Efb)) al- 
though a higher N c /N combined with a slightly lower E a 
would suffice to explain some of the data. We cannot there- 
fore exclude this mechanism but find that our present data 
suggest that the indirect chemisorption on SGs is more 
probable. 

There are several complications which we have ne- 
glected in the above treatment. One is the possib l e pres- 
ence of interstitial sites as suggested by IPulevI (^996) 
which can increase the efficiency of direct chemisorp- 
tion. Another perhaps is H2 formation through a reac- 
tion between two chemisorbed H atoms as suggested by 
ICazaux fc Tielensl l)2002[l2003t) . This mechanism would be 
efficient at high dust temperatures (>100K) and could be 
a possibility to explain the formation of H2 in Orion. 

In summary, we conclude that in order to explain the 
H2 formation efficiency in PDRs, the indirect chemisorp- 
tion mechanism upon small grains is the most promis- 
ing. This requires an activation barrier energy between a 
physisorbed H atom and a neighbouring chemisorbed H 
atom E a <1000 K (or <0.1 eV) and a fraction of occu- 
pied chemisorbed sites of around ten percent. This conclu- 
sion is consistent with our finding of a correlation between 
the H2 and PAH emission which suggests that Rf scales 
with the PAH abundance (see Sect.EJ). However, a better 
knowledge of the SGs properties (temperature, coverage 
of absorbed H atom, abundance ...) and of the mobilities 
of H atoms on realistic grain surfaces are critical. 

7. Conclusion 

The main aim of this study has been to provide estimates 
of the molecular hydrogen formation rate in a sample of 
nearby PDRs using results from both ISO and ground- 
based telescopes. The physical conditions in the PDR lay- 
ers from which H2 emission is observed (uh — 10 3 — 10 5 
cm~ 3 , T gas > 300 K) differ considerably from those in the 
diffuse clouds where one can observe the UV lines of H2 
(n H -100-1000 cm" 3 , T gas -50-100 K) and upon which 
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Fig. 5. Panel (a) : H2 formation rate as predicted for the indirect chemisorption mechanism (solid and dashed lines) 
and for the pure Langmuir-Hinshelwood mechanism (dotted lines) and as observed (crosses, see Table\^) as a function 
of the dust temperature. The dust temperatures for Oph W and S140 derived are similar (see Table\^j. The upper and 
lower panels show the rates predicted respectively for formation on small and big grains. The solid and dashed lines 
correspond respectively to E a equal to 1000 and 600 K. Panel (b) : H2 formation rate as predicted for the H2 formation 
by direct chemisorption (solid and dashed lines) and as observed (crosses) as a function of the gas temperature. The 
"observed" gas tem peratures are infer red from the distribution of low H2 pure rotational levels reported in Tabled for 
the Chamaeleon see lGrv et all \200A) . 



most estimates of the molecular hydrogen formation rate 
have been based. Thus the results from PDRs allow im- 
portant constraints to be placed upon the mechanisms for 
forming molecular h ydrogen on grain sur faces. We confirm 
the earlier result of iHabart et al.l l)2003|) that the H2 for- 
mation rate in regions of moderate excitation (\ <1000) 
such as Oph W, S140 and IC 63 is a factor of ^5 times 
larger than the standard rate estimated by Jura for dif- 
fuse clouds (and confirmed by recent FUSE data). On the 
other hand, towards regions of higher radiation field such 
as the Orion Bar and NGC 2023, we derive H2 forma- 
tion rates consistent with the standard value. Thus, the 
higher grain and gas temperatures in PDRs do not seem 
to impede the formation of H2 . 

We have attempted to interpret these results with sim- 
ple empirical models of the formation of H2 on grain 
surfaces. From these, we conclude that an "indirect 
chemisorption" m o del analogo us to that discussed by 
ICa.za.nx Tielens! l|2002l 12003ft is capable of explaining 
the data. This result requires an activation barrier en- 
ergy against the recombination of a physisorbed H atom 
and a neighbouring chemisorbcd H atom E a <0.1 eV. 
Another condition appears to be that one needs an ap- 
preciable fraction of surface sites occupied (few percent 
at least) with a binding energy of order 1 eV relative 
to the total number of surface binding sites (presum- 



ably mainly physisorbed with binding energies of order 
0.05 eV). Moreover, we suggest that small (size < 100 A) 
grains may be the most promising surface for forming H2 
in PDRs. There is in fact enough grain surface in small 
grains to allow the formation rate to be larger than the 
standard H2 formation rate and small grains spend most of 
their time at low temperature ijGuhathakurta fc DraineL 
H2 formation by indirect chemisorption upon small 
grains should be effective even in highly excited regions 
where large grains are quite warm. 

Our results show that formation of molecular hydro- 
gen in PDRs is likely to take place with a different mecha- 
nism than in the diffuse interstellar medium where the for- 
mation by physisorbed H atoms (Langmuir-Hinshelwood) 
probably dominates. In fact, in cold diffuse clouds the sur- 
face density of strongly bound H atoms should be low and 
consequently the formation by chemisorption would not 
be efficient. 

There are several fundamental uncertainties in our 
present estimates of the H2 formation rates in PDRs which 
future work should try to eliminate. One is due to the 
fact that we have used steady-state PDR models. This as- 
sumption can cause appreciable errors and it would be use- 
ful to calculate the expected H2 line intensities for models 
with lower H2 formation rates but where advection has 
been taken into account. However. IStorzer fc Hollenbachl 
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( 1998) have modelled PDR structure assuming an ion- 
ization front moving into the PDR and found that non- 
equilibrium effects are probably minor in objects similar 
to the Orion Bar. 

Other uncertainties are of an observational nature. It 
would be useful to have reliable estimates of the incli- 
nation and density distribution in the PDRs of interest 
in order to better constrain the m odels. There are indica- 

n i20on|) 



tions in the Orion Bar for instance ( Wal mslev et alJ 
that the column density is higher in the molecular layers 
of the PDR than in the region (discussed in this paper) 
where the H2 lines are formed. Density gradients perpen- 
dicular to the PDR photodissociation front clearly need 
to be taken into account when considering the spatial dis- 
tribution of the various H 2 lines. 

The theoretical models of H 2 formation discussed 
here are clearly very preliminary. More detailed models 
need to explicitly consider the degree of occupation of 
chemisorbed sites by H- atoms as well as the mobility of H 
atoms on various types of grain surfaces. Our results how- 
ever do suggest that more detailed consideration of H2 
formation on the surface of small grains would be worth- 
while. In this case, one should ideally follow the thermal 
fluctuations of these small particles and take into account 
the dependence on grain size in order to examine prop- 
erly their contribution to H2 formation. Clearly also our 
present estimates for processes such as tunneling are very 
crude. 

Additional information on the H2 formation process 
could probably be obtained using high quality data for H2 
in excited vibrational states which gives constraints on the 
excitation state of the newly formed H2. It also would be 
useful to obtain estimates of the H2 formation rate under 
a variety of different conditions. One such condition might 
be in the Magellanic clouds where the different metallic- 
ity, extinction curve, and radiation field potentially may 
influence both the available grain surface area and the 
efficiency of H2 formation. Another is in the thin h i gh ex- 
citation clouds such as that found by iMever et alJ l)200l[) 
towards the exciting star of NGC 2023. It is clear that 
given our uncertainty about grain compositions and size 
distributions in different ISM locations, one is forced to 
some extent to use the astrophysical data to guide our 
estimates of processes such as H2 formation. This is per- 
haps philosophically not as satisfactory as the traditional 
approach of employing experimentally determined or the- 
oretically calculated rates but it is likely that nature does 
not give us the choice. 
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Appendix : Probability for the diffusion and 
recombination of a physisorbed H atom with a 
chemisorbed H atom 



-rr (Eq. El see Sect. 



To calculate the probability f = - 

15.21 . we estimate the timescale r p for diffusion over the grain 
surface followed by recombination and the timescale r ev for 
evaporation of the physisorbed H. The latter time scale is given 
by 



r e „ = u Q exp 



( E d 



(8) 



with Ed the desorption energy of a physisorbed H and vq the 
vibrational frequency of H in a physisorbed site typically of the 
order of 10 12 s" 1 . t p can be given by the sum of the time to find 
a neighboring site to a chemisorption site bearing a H atom, 
r m and of the time to recombine and form H2, Tree- We assume 
that diffusion fro m one phys i sorbe d site to another occurs by 
thermal hopping jKatz et allll999h . The time to hop from one 
physisorbed site to the next is thus given by Uq 1 exp ( kr/ t ) 
with E p the activation barrier energy for physisorbed H atom 
diffusion. The mobility is a random walk and considering that 
there are 4 neighboring sites for each chemisorbed site, we find 
that one needs approximately (N/AN C ) 2 steps to be adjacent 
to a filled chemisorbed site. Thus: 



N 



x v exp 



f E p 

\ kTdus 



(9) 



The physisorbed H atom must cross the activation bar- 
rier to recombine with the nearest chemisorbed H atom 
by either (i) thermal diffusion with a probability f t h = 
exp (— fe7 f a - ) or by (ii) tunnelling with a probability f tU n = 
exp(-2-^£ (2m H E a ) - 5 ) jMessiahl IT97I with Ax the width 
of the barrier. The recombination time scale can be approxi- 
mated by 



N 
W c 



1 



{fth + ftun) 



(10) 



At high dust temperatures (above ~40 K) thermal diffusion 
dominates whereas at low temperature (below ~30 K) tunnel- 
ing is more important. 

There are many caveats to the above procedure. 
In particular, we note that considering the probability 
exp (— 2 (2mHE a ) ' 5 ) for tunneling transmission through 
a single barrier is incorrect for relevant astrophysical surfaces 
which are not regular and need therefore to be modeled by 
a set of barriers as many as there are atoms. The tunneling 
probability we adopt is thus an upper limit. However, for the 
purpose of comparison with observationally determined rates, 
our approach allows us to see which of the H2 formation pro- 
cesses considered in this study will dominate. 
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